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ON SIMILARITY RULES FOR TRANSONIC FLOWS 

By Cabl Kaplan 


SUMMARY 

A method used by Tsien to derive similarity rules for hyper- 
sonic flows is utilized to derive Yon Karmdn’s similarity rules 
for transonic flows. A slight generalization is introduced by 
the inclusion of y, the ratio of specific heats, as a parameter. 
At the lower limit of the transonic region of flow the theory 
yields a formula for the critical stream Mach numbers of a 
given family of symmetrical profiles. It is further shown, that 
this formula can also be obtained by means of the Prandtl- 
Olauert small-perturbation method. Investigation of the be- 
havior of the similarity parameter in the region where the 
thickness coefficient approaches zero and the critical stream 
Mach number approaches unity shows that it possesses a 
limiting value characteristic of the prescribed family of shapes. 

INTRODUCTION 

Tlie rigorous solution of the subsonic flow of a compressible 
fluid past a prescribed closed body thus far has proved to 
be of insurmountable difficulty. As a consequence of this 
difficulty the emphasis has been placed on the establishment 
of a correspondence between the flow past a given body in an 
incompressible fluid and the same body in a compressible 
fluid. Among the best known results of this mode of attack 
are the Prandtl-GIauert rule and the Yon Kfirm&n-Tsien 
velocity or pressure correction factor — both based on some 
form of linearization of the fundamental nonlinear flow equa- 
tions. None of the methods based on the linearization of the 
flow equations, however, can yield correct results in the 
transonic range where the flow is partly subsonic and partly 
supersonic. For this case a certain amount of the feature of 
nonlinearity of the flow equations must be retained in order 
to obtain useful and nontrivial results. In the present paper 
a detailed derivation is given of the transonic similarity 
rules recently given by Yon Kdrm&n (reference 1 ). 


FUNDAMENTAL EQUATIONS 

In plane steady flow the equation governing the flow of a 
nonviscous compressible fluid can be written in the form 


li-” (5+5)+ <«’-'•) £r 0 « 


In the derivation of this equation the pressure is assumed to 
be a function of the density only. If, further, the motion is 
irrotational, then 


dy dr 


( 2 ) 


Here u and v are, respectively, the component velocities along 
the Cartesian x- and y-axes and c is the local velocity of 
sound given by 

c *=c 0 2-^=±(u , + s ! ) (3) 


where c 0 is the velocity of sound at a stagnation point 
u=0, v=0 and 7 is the ratio of specific heats at constant 
pressure and constant volume. 

Equation ( 1 ) is far too complicated to afford an insight into 
the properties of potential Sow in the neighborhood of 
Mach number unity. The discussion is, therefore, restricted 
to the flow past a thin profile. Thus, at first 0 is assumed 
to be small in comparison with the sound velocity c 0 . Equa- 
tion ( 1 ) is then simplified to 


7+1 u 3 
2 c? du do 
7—1 u 3 dx~^dy~ 

2 Co* 


(4) 


By the introduction of the sound velocity for which the 
local fluid and sound velocities are equal 



(5) 


and of the maximum possible fluid velocity 



2 

7—1 


Co 


equation (4) becomes 

c* 2 C>M , dr 
, u* dx^~dy 
1 a 3 

Umax 


( 6 ) 


(7) 


It is desirable to. simplify this equation still further but yet 
to retain those features which yield nontrivial and useful 
results. Thus, if it is assumed that u is of the order of the 

critical velocity c* and if only terms up to the order 1 — ^ 
are retained, then equation (7) can be mitten as 

-(7+l)(«-e*)g+c*g=0 (8) 

Thus far the irrotationality condition, equation ( 2 ), has 
not been used. If, now, the undisturbed stream past a 
slender body is of velocity U slightly different from the 
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velocity of sound and in the direction of the positive z-axis, 
then, according to equation (2) and the assumptions leading 
to equation (8), a velocity potential $ can be introduced with 


$ — c*x+ (1 — M*) <p 


(9) 


TJ ' " . 

where M *=— * and (1— M*)<p is the disturbance-velocity 
c 

potential. 

Tlien 


sf 1 


( 10 ) 


and equation (8) becomes 




( 11 ) 


Equation (11) is a nonlinear simplified form of . the fundamen- 
tal differential equation (1) and has been treated recently 
by Von Kdrm&n in connection with similarity rules in two- 
dimensional transonic flow (reference 1). Equation (11), 
when expressed in hodograph variably, is of the type treated 
by the Italian mathematician F. Tricomi some years ago 
(reference 2). 

DERIVATION OF SIMILARITY CONDITIONS 

Recently, Tsien (reference '3} derived similarity rules for 
hypersonic flows where the fluid velocity is much larger than 
the velocity of sound. In the present paper the same proce- 
dure is employed to derive similarity rules for transonic flow's 
where the fluid velocity is very nearly that of sound. 

According to the assumptions leading to equation (11), it 
is implied that the solid body is thin and possesses no stag- 
nation points and that the velocity of the fluid is everywhere 
in the neighborhood of the local velocity of sound. Now, 
suppose the profile of the obstacle to be symmetrical with 
respect to both the z- and y-axes and to possess cusps at both 
the leading and the trailing edges. Such profiles with uni- 
form flow in the direction of the long axis-of-9ymmetry x 
fulfill the assumptions leading to equation (11). The flows 
past these profiles are said to be similar if the equation of 
motion (11) and the boundary conditions can be expressed 
in nondimensional variables in such a way that only a single 
constant factor is involved. Thus, if 2a is the chord and 26 
is the maximum thickness of the body, then the following 
nondimensional variables are introduced: 


x=a£ ^ 

y=a(y+l) m t n i}j 


( 12 ) 


where t=— and m and n are exponents yet to be determined. 
0/ 


It is clear that the nondimensional quantity involved is the 
thickness coefficient t since this quantity determines the mag- 
nitude of the disturbance velocities. The exponents m and n 
are to be determined in such a w r ay that the same constant 
factor appears in both the equation of motion and the 
boundary conditions. 

The appropriate nondimensional form for the velocity 
potential <p is 

<p=ac*f(Z,r)) (13) 

By substitution from equations (12) and (13), the equation 
of motion (11) becomes 

- (Y+ 1) J » +l (1 - M*) 0 (14) 

The boundary conditions at infinity require that the flow 
velocity be U. Hence, from the first of equations (10), 


<V = * 

bar 

^=0 

by 


(at oo ) 


The boundary condition at the surface of the slender body is 
c*g (£) t (at y= 0, - a Z x<a) 

where g describes the distribution of ’slope along the 
surface of the body. 

By means of equations (12) and (13), the boundary con- 
ditions can be weritten as: 


I — 1 


V (at co) 


(7+l)-’»(l-M*)r- 1 ^=y(© (atu=0,-l£«l) 

(15) 

A comparison of the differential equation (14) and the bound- 
ary conditions, equations (15), shows that a single parameter 
is involved if 

2n=^-(n+l) or n=~ i 


that is, 


2m+l = ^m or m=— ^ 


( 7 +l) 1/8 (l-M*)r 2 '*= 2 E: 
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The undisturbed-stream Mach number M m —— can be 

introduced in the following way: 

The general relation between M* and M„ is 


... ( 4 ^ Y 


(16) 


or 


i/i 


(1+^ (l-flf*)[2-(l-M*)]j 
Then, if powers of l—M* higher than the first are neglected, 


Therefore 


1— M*— (1— M a )-\r 


(l-Jl/„)[(T+l)fl- 2/8 =Z 


(17) 


The results obtained thus far are such that by means of the 
substitution equations 


x—a.% 

y=a[(y+l)t]~ ul v 
¥>=ffc*/(l, v) 

#=(l-fif„)[(7+l)r 2 ' !! 


(18) 


the differential equation for /(£, 17 ) and the boundary con- 
ditions become: 




(19) 


and 


^=-1 

?-0 




(at 00 ) 


(at v=0,— l^l^l) 


( 20 ) 


The meaning of the similarity rule implied in the defi- 
nition of the parameter K is the following: 

If a series of bodies having the same distribution function 
< 7 (f) for the slope but different thickness ratios t are placed 
in flows of different undisturbed-stream Mach numbers 
M„ and different values of 7 , such that the parameter 

l—M 

K— [(y-f,' ;[ )fp /3 remains constant, then the flow patterns are 

similar in the sense that the same function /(£, 17 ) describes 
the flows. 


RESULTS DERIVED FROM THE SIMILARITY RULE 

PRESSURE COEFFICIENT 

In the case of a uniform flow past a fixed boundary, the 


pressure coefficient is defined as 


Op, n a 


P-Jpo 




where p a and p„ are, respectively, the pressure and density 
in the undisturbed stream and the static pressure p in the 
fluid is given by 


P=P~ [ 


1 


(t-i) r ™ 


Then 

a 


Pi M a 




1-4 


7—1 




( 22 ) 


By means of equations (10) and (18), if powers of 1 — M a 
higher than the first are neglected, the following result is 
obtained: 

C p ,ir a — (1 M a ) 
or 


t 1 ' 3 


C’PliTa 1) !/3 ■?*(£> V i 


(23) 


where P(£, 17 ; K) depends on the form of the solution /(£, 17 ) 
for the particular family of profiles treated. 


LEFT COEFFICIENT 


The lift l of the body is given by 
(=(j) (p),_ oda: 


By a similar procedure, as in tlie derivation of equation (23), 
the lift coefficient is given by 


c t = 


l 


v V3 


| P „Z7 2 (2a) ( ' y+1 )' /3 


L{K) 


(24) 


where 


L{K) = J* ^ P(£, 0 ',K)d% 


and where for an extremely thin straight-line profile the 
thickness coefficient t has been replaced by the angle of 
attack a. 

DRAG COEFFICIENT 

The pressure drag d of the body is given by the following 
expression: 

d=2 f-a AX 
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Hence the drag coefficient is given by 


Hence equation (29) becomes 


*- T—^ <25) 

\ P .v'm l7 ) 

where 

D(K)=2§\{i)P{t,0\K)dZ 

ADDITIONAL CONSIDERATIONS 

The results derived in the present paper apply to two- 
dimensional near-sonic flows past thin shapes. .Such flows 
have been calculated for a family of symmetrical shapes 
with cusped leading and trailing edges (reference 4) and for 
a family of elliptic cylinders (reference 5). These calcula- 
tions are valid at least up to the critical stream Mach num- 
ber M cr . The critical Mach number may be considered 
from two points of view. First, it may be considered, to 
denote the lower limit of a mixed subsonic-supersonic flow, 
that is, where the imbedded supersonic region is simply the 
point of maximum fluid velocity at the surface of the solid. 
From this point of view, according to equation (17), a 
critical value of the similarity parameter K can be defined. 
Thus 

Kcr= [(y+i)t] vi ' ' C26) 


MJ 


a -M„r n 




( 30 ) 


As seen from equation (29), 1 —M C r is of the first order in the 
small perturbation u'/U and, accordingly, equation (30) 
includes terms of higher order. This fact can be seen by 
rewriting equation (30) in the following form: 


7+1 


(1— M c , 2 ) 8/2 




F 


7—1 

7+1 


(1 -MJ) 


' J >,0 


Then, to within the lowest order in the small perturbation, 


or, approximately, 

(31) 


It is quite easy to show the connection between equation (31) 
and equation (27). Thus, for an incompressible fluids the 
pressure coefficient is given by 

C?p,o— 1 - (32) 


This equation can also be •written in a form that yields the 
critical stream Mach number for a given family of shapes; 
that is, 

M cr =l-K cr [(y-+l)t] i '* (27) 

Second, the critical-stream Mach number may be consid- 
ered to denote the upper limit of the purely subsonic range 
of speeds. This point of view suggests a derivation of equa- 
tion (27) by means of the Prandtl-Glauert small-perturbation 
method (reference 6). The procedure is as follows: 

The relation between the local and the undisturbed-stream 
Mach number, within the approximation of the small- 
disturbance theory, is given by 

M s =M 0> s [l+2^(l+^M„ 2 )] (28) 

where u' is the disturbance velocity and 17 is the undisturbed- 
stream velocity. 

fjJ 

By definition, M— 1 for = M cr and, §ince C f , = — 2 jj - 

approximately, equation (28) becomes 

1 =MJ [l - C v , M cr (l+^p w)] (29) 

The relation between C p , u „ and the pressure coefficient C v , 0 
of the incompressible fluid is given by the Prandtl-Glauert 
rule 

n o 

yi -m„ 2 ' ' ' 


where q is the magnitude of the fluid velocity at any point in 
the field of flow. 

In tbecase of the family of shapes of reference 4, the maxi- 
mum velocity at the surface is given by 


#-1+1*+ ■ ■ 


Hence 


0 Pi a == — 3t+ . . . 
and equation (31) becomes 


M. 


- l ~Q 


32; 


1/* 


[(7+1 )tY* 


(33) 


In the case of the family of elliptic cylinders of refer- 
ence 5, the maximum velocity at the surface is given by 


TT 1+t 


Hence 


C Pt o= — 2<+ . 
and equation (31) becomes 


M cr =l-|[(7+l)tr 


(34) 


An examination of equations (33) and (34) suggests that 
K„ possesses a limiting value; that is, 


lim [( * + ^p 

M.r-il 


( 35 ) 
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It is noteworthy that the value of (K„) lim depends on the 
family of shapes, although in the limit f-*0 the profile in 
every case is a straight-line segment. The numerical values 
of (Kcr) l{m shown in equations (33) and (34) represent, how- 
ever, only the effect of the Prandtl-Glauert or first-order 
term in the power-series development in t of the maximum 
velocity at the solid surface. It is rather surprising that the 
higher-order terms involving the higher powers of t also 
contribute to the value of (K„) l(m . The procedure is simply 
to replace the maximum velocity 'by the critical speed; that is, 



where a u a t , a t , . . . depend on the given f amil y of shapes and 
involve only y and M„. 

Then as suggested by equation (26), when 1 —M„ 2 is re- 
placed by 2K CT [(y+l)t\ v \ equation (36) is identically satis- 
fied for t— >0, with each additional higher-order term con- 
tributing to the value of K„. In this fashion there are 
obtained, successively, linear, quadratic, cubic, and higher- 
order equations for the determination of (K„) Uni . 

The foregoing considerations have been applied to the 
family of symmetrical shapes of reference 4 and to the family 
of elliptic cylinders of reference 5. The values of K„ for 

TABLE I . 

VALUES OF K„ FOR THE FAMILY OF SYMMETRICAL 
SHAPES OF REFERENCE 4 
[For air, 7=1.4; lor Freon-12, 7=1.136] 


f 

Afcr 


Approximation. 

Approximation' 

First 

Second 

Third 

First 

Second 

Third 

Air 

0 

.002 

.005 

.010 

.020 

.040 

.060 

.080 

.100 

m TyW 

|B 

1.000 

.980 

.963 

.942 

.909 

.859 

.818 

.784 

.753 

wWTmm 

(1635 

.647 

.638 

.632 

.624 

.800 

.393 

.581 

.568 

a 717 

.710 

.705 

.700 

.689 

.675 

.664 

.650 

.639 

0.738 

.732 

.726 

.719 

.709 

.696 

.685 

.678 

.065 

Freon-12 

0 

.002 

.005 

.010 

.020 

.040 

.060 

.080 

.100 

L 000 
.083 
.960 
.951 
.924 
.881 
.848 
.819 
.794 

1.000 

.981 

.966 

.946 

.915 

.887 

.830 

.797 

.709 

m nju 

BE 

IWjrr’iM 

0.717 
.713 
.709 
.704 
. 694 
.682 
.670 
.65S 
.647 

0.738 

.735 

.730 

.725 

.716 

.706 

.60S 

.691 

.682 


values of t different from zero were obtained from equa- 
tion (26) with the required values of M u determined by means 
of the theoretical results of references 4 and 5. The limiting 
values of K a were obtained from equations corresponding to 
equation (36). Tables I and II show the results of these 
calculations for air and for Freon-12. Note the approach 
to the limiting values of K„ as t approaches zero and M„ 
approaches unity. The first column of values of K„ shows 
the effect of the first-order or Prandtl-Glauert term; the 
second and third columns show, respectively, the effects of 
the second-order and third-order terms in the thielmftga 
coefficient t. The successive values of K cr for the various 
values of the thickness coefficient t indicate good convergence. 

Langley Memorial Aeronautical Laboratory, 
National Advisory Committee for Aeronautics, 
Langley Field, Va., December 9, 1947. 
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TABLE II 

VALUES OF K„ FOR THE FAMILY OF ELLIPTIC CYLINDERS 
OF REFERENCE 5 


[For air, 7-1.4; for Freon-12, 7= 1.136] 


C 

M, r 

AL, 

Approximation 

Approximation 

ESI 

Second 

Third 

First 

Second 

Third 

Air 

0 

.002 

.005 

.010 

.020 

.040 

.060 

.080 

.100 

LOGO 

.986 

.974 

.960 

.936 

.901 

.873 

.848 

.827 

1. 000 
.985 
.972 
.957 
.932 
.894 
.863 
.837 
.813 

1. 000 
.984 
.071 
.955 
.9®} 
.892 
.861 
.833 
.809 

a 500 
.494 

.490 

.487 

.481 

.470 

.462 

.455 

.417 

a 637 
.530 
.626 
.522 
.515 
.505 
.497 
.490 
.484 

0.556 

.650 

.514 

.537 

.526 

.515 

.507 

.500 

.494 

Freon-12 

0 

.002 

.006 

.010 

.020 

.040 

.060 

.080 

.100 


LOCO 

.980 

.974 

.959 

.936 

.900 

.872 

.847 

.826 

B ml 

B 

B :-m» 

& 

a 537 
.534 
.532 
.528 
.520 
.509 
.504 
.496 
.488 

0.556 

.553 

.550 

.543 

.534 

.522 

.513 

.606 

.498 












































